Fractional Quantum Hall Effect of Hard-Core Bosons in Topological Flat Bands 
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Recent proposals of topological flat band (TFB) models have provided a new route to realize the 
fractional quantum Hall effect (FQHE) without Landau levels. We study hard-core bosons with 
short-range interactions in two representative TFB models, one of which is the well known Haldane 
model (but with different parameters). We demonstrate that FQHE states emerge with signatures 
of even number of quasi-degenerate ground states on a torus and a robust spectrum gap separating 
these states from higher energy spectrum. We also establish quantum phase diagrams for the filling 
factor 1/2 and illustrate quantum phase transitions to other competing symmetry-breaking phases. 
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Introduction. — The fractional quantum Hall effect 
(FQHE), one of the most fascinating discoveries in two- 
dimensional (2D) electron gas, has set up a paradigm 
to explore new topological phases in other strongly cor- 
related systems. As commonly believed, the FQHE re- 
quires two basic ingredients: single-particle states with 
nontrivial topology, and quenching of the kinetic en- 
ergy compared to interaction energy scale. However, de- 
spite of the seemingly universal theoretical concepts, the 
FQHE has only been found in 2D systems under a strong 
perpendicular magnetic field, i.e., in which particles move 
in Landau levels (LLs). In rotating Bose-Einstein con- 
densate [lj and optical lattice systems 0, 0] , researchers 
have been interested in generating an artificial uniform 
magnetic field, thus the bosonic FQHE states are ex- 
pected, but still due to the existence of LLs. 

Haldane's honeycomb lattice model [4J and other sim- 
ilar lattice models 0, Q have two nontrivial topological 
bands characterized by ± 1 Chern numbers @, Q j demon- 
strating the integer quantum Hall effect without LLs. 
However, these single-particle bands are still highly dis- 
persive, thus it is unlikely to realize the FQHE in such 
systems. Recently, proposals of topological flat bands 
(TFB) [§-[ll| shed new lights to overcome this long- 
standing and hard problem. These TFB models belong 
to the same topological class as the Haldane model and 
are distinct from other fiat bands with a zero Chern num- 



ber [12j . The energy dispersion in such a TFB is substan- 
tially reduced by tuning short-range hopping parameters. 
In particular, based on the mechanism of quadratic band 
touching jj| [ll| , a series of TFB models have been ex- 
plicitly constructed with a flatness ratio (the ratio of the 
band gap over bandwidth) reaching a high value between 
20 - 50. 

Here, we focus on the possible bosonic FQHE in TFB 
models filled with interacting hard-core bosons, since the 
TFB will be more likely realized in optical lattices by ma- 



nipulating bosonic cold atoms 13|-15||. Although TFB 
models possess both ingredients to realize the FQHE, 
quantum phases in such systems are determined by some 
competing effects, different from a LL problem. The 
main effects are: i) the lattice effect and the residual ki- 
netic energy since a TFB is not strictly flat; ii) the Berry 
curvature of a TFB has substantial momentum depen- 
dence representing a non-uniform magnetic field effect 
in momentum space; iii) lattice symmetry breaking may 
lead to other conventional ordered states for hard-core 
bosons. Thus it is highly desirable to pursue a system- 
atic numerical study of such interacting TFB models. 

In this letter, we present the exact diagonalization 
(ED) calculations of two representative TFB models with 
the nearest-neighbor (NN) and the next-nearest-neighbor 
(NNN) repulsions V\ and Vi. We find convincing nu- 
merical evidences of both the 1/2 and the 1/4 bosonic 
FQHE phases which are characterized by the formation 
of quasi-degenerate ground-state manifold (GSM) with 
even number of states. The GSM carries a unit total 



Chern number [lfjfl , which is a robust property of the sys- 
tem protected by a finite energy spectrum gap. We also 
determine phase diagrams for our systems and illustrate 
the quantum phase transitions based on the calculations 
of the density structure factors and the fidelity [l7[ of the 
ground state (GS) wavefunction. 

Model Hamiltonians. — The first model is the Haldane 
model [J| on the honeycomb (HC) lattice filled with in- 
teracting hard-core bosons: 
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FIG. 1: (color online), (a) The Haldane model on the honey- 
comb lattice and (b) the checkerboard model. The arrow di- 
rections present the signs of the phases ±<j> in the NNN or NN 
hopping terms. For the checkerboard lattice, the NNN hop- 
ping amplitudes are t' (—t') along the solid (dotted) lines. The 
NNNN hoppings are represented by the dashed curves in both 
models, (c)-(d) Intensity plots of spectrum gaps in the V1-V2 
phase space at v = 1/2 for (c) 24-site honeycomb lattice; (d) 
24-site checkerboard lattice. FQHE, SF, SS1/SS2, and Solid 
label estimated phase regions inferred from the spectrum-gap 
plots and other information (see text). 



where b\ creates a hard-core boson at site r, n r is the bo- 
son number operator, (...), ((•••)) and (((. . . ))) denote 
the NN, the NNN and the next-next-nearest-neighbor 
(NNNN) pairs of sites, respectively [Fig. QJa)]. We call 
this model as a Haldane-Bose-Hubbard (HBH) model. 

The second model is a variant version of the HBH 
model on the 2D checkerboard (CB) lattice 0, [HL EM : 
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Topological flat bands. — On the honeycomb lattice, if 
we restrict the model with only NN and NNN hoppings, 
the best TFB has a flatness ratio 7 [10]. By allowing 
the NNNN hoppings, we numerically found several flat- 
ter bands with nonzero Chern numbers, e.g. a flatness 
ratio about 50 for the set of parameters, which will be 
used here: t = 1, t' = 0.60, t" = -0.58 and 4> = 0.4tt. Us- 
ing these parameters, the TFB is gaped from quadratic 
touching [19( by breaking the time reversal symmetry but 
preserving other lattice symmetries [llj . 

On the checkerboard lattice, we adopt the parameters 
of Ref. ll| with an additional minus sign (to make the 
TFB as the lower energy band), t = -1, t' = 1/(2 + y/2), 



t" = -1/(2 + 2\/2) and = tt/4, which leads to a TFB 
with the flatness ratio about 30. 

The v = 1/2 phase diagrams. — In the ED study, we 
consider a finite system of Ni x N% unit cells (total num- 
ber of sites N s — 2 x Ni x N2) with periodic boundary 
conditions (PBC), denoting the number of bosons as N b , 
and the filling factor of the TFB is thus v = N b /{N 1 N 2 ). 
In both models, the amplitude of NN hopping |f| is set 
as the unit of energy. We first glance at the spectrum 
gaps of the two 24-site (2x4x3) lattices at the filling 
v = 1/2 as shown in Fig. QJc) andQJd). E\, E 2 and E 3 
denote the energies of the three lowest eigenstates. For 
the v = 1/2 FQHE phase at the left bottom corners in the 
V1-V2 space, there is a GSM with two quasi-degenerate 
lowest eigenstates, and the GSM is separated from higher 
eigenstates by a finite spectrum gap E 3 — E 2 (^> E 2 —Ei). 
The other rough phase regions for the possible superfluid 
(SF), the supersolids (SS1/SS2) and the solid will be dis- 
cussed later. We have also obtained numerical results 
from larger lattice sizes of 32 (2 x 4 x 4), 36 (2 x 6 x 3) 
and 40 (2 x 4 x 5) sites, and have confirmed both phase 
diagrams are qualitatively correct. 




FIG. 2: (color online). 1/2-FQHE spectrum gaps versus 1/N S : 
(a)-(b) honeycomb lattice; (c)-(d) checkerboard lattice. 

Low energy spectrum and robust spectrum gap. — We 
denote the momentum vector q = (2nki/Ni,2irk 2 /N 2 ) 
with (k\,k 2 ) as integer quantum numbers. The GSM 
is defined as a set of lowest states separated from other 
excited states by a finite spectrum gap. If (hi, k 2 ) is the 
momentum sector for one of the state in the GSM, we 
find that the other state should be obtained in the sector 
(fci + N b , k 2 + N b ) [module (N h N 2 )} . For N s = 24, 36 
and 40, the two states within the GSM of a v = 1/2 
FQHE phase are indeed in different momentum sectors: 
(0, 0) and (2, 0) for N s = 24 and 7V S = 40, while (0, 0) 
and (3,0) for N s = 36. For N s = 32, both N b /Nx and 
N b /N 2 are integers, thus both states within the GSM are 
in the (0, 0) sector. 

Now we check whether the spectrum gap E3 — E 2 re- 
mains in the thermodynamic limit. As shown in Fig. [5J 
when N s increases, the spectrum gap E3 — E 2 does not 
decrease, which extrapolates to a finite value at large 
N s limit. Interestingly, the spectrum gap E3 — E 2 is 
already quite large (E3 — E 2 3> E 2 — E±) for hard-core 
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boson system without additional interactions (Vi=V2=0) 
[in Figs. (Ha) and HJc)] demonstrating the robust 1/2 
FQHE. The spectrum gap can be slightly enhanced with 
a small Vi and/or V 2 [Figs.^b) and^d)]. 

By comparing the spectrum gap E3 — E2 between two 
lattices for the Vi = V% = 0.0 cases in Figs.dja) and[5Jc), 
we notice that the gap E^~E2 in the honeycomb lattice is 
obviously larger, which might be due to its larger flatness 
ratio. After studies on a few more cases with smaller 
flatness ratios 7 Qlj and 30 (with t' = 0.40, t" = -0.33 
and 4> — 0.5ir) on the honeycomb lattice, we conclude 
that the flatter the TFB is, the larger the spectrum gap 
E3—E2 can be, indicating a more robust FQHE, although 
the global structure of the phase diagram does not change 
much. 

(a) HC, N s =24, V,=1 .0, V 2 =0.4 (b) HC, N s =32, V,=0.6, V 2 =0.2 
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FIG. 3: (color online), (a)-(b) Low energy spectra versus 81 
at a fixed 82 — for the honeycomb lattices at v = 1/2. (c) 
F(6i, 82) of a GSM of the 32-site honeycomb lattice. 



Berry curvature and Chern number. — Introducing two 
boundary phases 6\ and 62 as the generalized bound- 
ary conditions in both directions, the Chern num- 
ber (is also the Berry phase in units of 2-k) of a 
many-body state is given by an integral in the bound- 
ary phase space H Gil: c = h 1 1 dei9 2 F(9-L,9 2 ), 
where the Berry curvature is given by F(6i,02) = 



Im 



((«■ 



<9M> 

,)!)., 



For each GSM of 1/2- 

FQHE phase with N s = 24, 36, 40, the two states 
are found to evolve into each other with level crossing 
when tuning the boundary phases [Fig. |3Ja)]. While for 
N s =32, with both states of the GSM in the (0,0) sec- 
tor, each state evolves into itself when tuning the bound- 
ary phases, and avoided level crossings appear instead 
[Fig. [3][b)] due to nonzero coupling between same mo- 
mentum states. The GSM in the FQHE phase also has 
rather smooth Berry curvature [Fig. 03c)] and shares a 
total Chern number C = 1. 

SF stiffness and structure factors. — The 1/2 FQHE 
phase on the honeycomb lattice is also distinguished from 
the other phases by the featureless intra-sublattice (AA) 
structure factor 5 AA (q) [Fig. Ufa)]. The solid phase at a 
larger V 2 is characterized by a ridge with qi + q 2 = 27r in 
the inter-sublattice (AB) structure factor 5 AB (q) [Fig. 0] 
(c)] and an almost vanishing psF [Fig. life)]. The SF 
phase at a smaller V2 has the finite psf [Fig. IUe)] but 
with a weaker ridge in 5 AB (q) [Fig. @] (b)]. At a fixed 




FIG. 4: (color online). 32-site honeycomb lattice at v = 1/2. 
(a) S AA (q) of the FQHE phase; (b) S AB (q) of the SF phase; 
(c) S AB (q) of the solid phase, (d) Excited energy E„ — Ei in 
various sectors (with d-fold degeneracy) versus V2 at V\ — 4.0. 
(e) psf, S ab (ty,tt), GS overlap |(*(^ 2 )|*(10))|, and fidelity 
susceptibility xf versus V2, at a fixed V\ = 4.0. (f) Illustration 
of boson occupancy in the solid phase. 



Vi = 4.0 while tuning V2, a transition from the FQHE 
to the SF phase occurs with the level crossing of E 2 and 
£3 around V2 — 1.0 [Fig. |4jd)] ; and a transition from the 
SF phase to the solid phase near V2 = 2.5 is indicated by 
a peak of the GS fidelity susceptibility [ijj in Fig. Bfe), 
where X f = 2[1 - |(*(V r 2 )|*(^ 2 + S)}\}/6 2 . 




V,(=V 2 ) 



FIG. 5: (color online). 32-site checkerboard lattice at v — 1/2. 
(a) 5" AA (q) of the FQHE phase, (b) 5* AB (q) of the SSI phase; 
(c) S AA (q) of the SS2 phase, (d) p SF , S AA (7r,7r), and X f 
versus Vi along the Vi = V2 line, (e)-(f) Illustration of boson 
occupancy in the SSI and SS2 phases, respectively. 

Similarly, the 1/2 FQHE phase on the checkerboard 
lattice also differs from the other phases by feature- 
less 5(q) [Fig. Eta)]. While both SSI and SS2 phases 
are characterized by either the q = (tt/2,tt/2) peak of 
5 AB (q) [Fig. or the q = (tt.tt) peak of S* AA (q) 

[Fig. Etc)]. Along the V x = V 2 line while V x {= V 2 ) be- 
ing tuned, a transition can be inferred from the FQHE 
phase to the SS2 phase around V\(= V 2 ) = 1.0 by a 
sharp peak in the GS fidelity susceptibility [Fig. EJd)]. 
We emphasize that the firm establishment of the super- 
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solid and solid phases needs scaling of both psF and 
5(q) for systems with larger sizes, e.g. 2x6x6 and 
2x8x8, which are compatible with the ordering patterns 
[Figs. Ilff), [5je), [5jf)] but are far beyond the capability 
of the present ED method. 
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FIG. 6: (color online). Fhe 1/4-FQHE of the 40-site checker- 
board model: (a) spectrum gaps E n — E\ versus k\N2 + fe; 
(b) low energy spectra versus 9\ at a fixed 82 = 0. 

The v = 1/4 FQHE. —We searched for the v = 1/4 
FQHE for both lattice models with the system sizes of 
N s = 24, 32, 40 and 48. In contrast to the v = 1/2 
case, the onset of the 1/4-FQHE needs finite values of 
either V\ or Examples of the 40-site checkerboard 
system are chosen to demonstrate the basic properties of 
the 1/4-FQHE. For each set of V x and V 2 in Fig. EJa), 
there is clearly a GSM with four states. For each GSM, 
the four states evolve into each other when tuning the 
boundary phases [Fig.^b)], and all the four states share 
a total Chern number 1. These are concrete evidences of 
the 1/4-FQHE in some small parameter regions, however, 
these regions depend on the lattice sizes more sensitively 
than those of the 1/2-FQHE. We conjecture that a finite 
NNNN repulsion V3 may be necessary to get a large and 
stable parameter space of the 1/4-FQHE, which will be 
addressed in a future work. 

Summary and discussion. — We consider hard-core 
bosons in two representative TFB models with NN and 
NNN repulsions. We find convincing numerical evi- 
dences of both the 1/2 and the 1/4 bosonic FQHE phases 
which are characterized by distinctive finite spectrum 
gap; quasi-degenerate states in a GSM which can evolve 
into each other upon varying boundary phases; smooth 
Berry curvature and topologically invariant unit total 
Chern number for the GSM. For both lattices, the 1/2- 
FQHE phase is found to occupy a significant space of 
phase diagrams, in addition to other conventional or- 
dered phases. Interestingly, such a 1/2-FQHE is very 
stable (large spectrum gap) for hard-core bosons even 
without additional interactions (V\ = V2 = 0), which 
makes it easier to be realized by cold atoms in optical 



lattices. 

Our work here focuses on hard-core bosons in TFB 
models, while a recent parallel work found the 1/3 and 
1/5 FQHE of interacting fermions on the checkerboard 
model [2^] (see also an updated version of Ref. [13). 
There is also an interesting proposal lately to find the 
bosonic FQHE (in terms of chiral spin states) in frus- 



trated kagome-lattice magnets [21 
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